Correlation energy of the one-dimensional Coulomb gas 
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We introduce a new paradigm for finite and infinite strict-one-dimensional uniform electron gases. 
In this model, n electrons are confined to a ring and interact via a bare Coulomb operator. In the 
high-density limit (small-r s , where r 3 is the Seitz radius), we find that the reduced correlation energy 
is e c (r a ,n) = e' 2 '(n) + 0(r s ), and we report explicit expressions for e' 2 '(n). In the thermodynamic 
(large-n) limit of this, we show that e c (r s ) = — 7r 2 /360 + 0(r s ). In the low-density (large-r s ) limit, 

the system forms a Wigner crystal and we find that e c (r s ) = — [ln(v / 27r) — 3/4] rj 1 +0.359933 r^ s ^ 2 + 
0(r^ 2 ). Using these results, we propose a correlation functional that interpolates between the high- 
and low-density limits. The accuracy of the functional for intermediate densities is established 
by comparison with diffusion Monte Carlo results. Application to a non-uniform system is also 
reported. 



Introduction. — The usual paradigm for modeling uni- 
form electron gases (UEG) in one dimension (ID) is 
the well-established Luttinger model [l[ , which is exactly 
solvable by a Bogoliubov transformation technique [2j 
and is particularly useful to study the low-energy spec- 
trum of the ID UEG This model has been moti- 
vated by the search for a replacement for the Fermi liq- 
uid theory in ID which fails because of some divergences 
and makes conventional fcrmion many-body perturbation 
theory impracticable. The Luttinger model is shown to 
provide a general method for resumming all of the di- 
vergences [4j . Although features associated with the Lut- 
tinger model have been observed in various systems, such 
as semiconductor quantum wires, ultracold atoms con- 
fined in elongated traps, GaAs quantum wells, carbon 
nanotubes, and many others Q , the model is not strictly 
applicable to these systems because, whereas the stan- 
dard Luttinger model assumes a short-ranged interaction, 
electrons actually interact via the long-range Coulomb 
force. 

In the literature, the bare Coulomb operator 1/x 
(where x is the interparticle distance) is usually avoided 
because of its divergence at x = and the intractability of 
its Fourier transform in ID. Instead, most studies of the 
ID UEG adopt a quasi-lD description by adding a trans- 
verse harmonic potential 043 ot use a potential of the 



form 1 



jjfl , where \i is a parameter that eliminates 



the singularity at x = while retaining the long-range 
Coulomb tail 0-@|- However, the introduction of the pa- 
rameter fj, is undesirable, for it modifies the physics of 
the system, especially in the high-density regime where 
neighboring electrons experience only the constant part 
of the potential. It is also unnecessary because, if the 
wave function vanishes when any two electrons touch, 
the Coulomb potential does not lead to energy diver- 
gence This allows us to apply the ID Bose-Fermi 



divergence of the Coulomb potential has the effect of 
mimicking the Pauli principle, which, for fermions, pro- 
hibits any two fermions from overlapping. This implies 
that, for strictly ID systems, the bosonic and fcrmionic 
ground states are degenerate and the system is "spin- 
blind" . Consequently, the paramagnetic and ferromag- 
netic states are degenerate and we will consider only the 
latter Q. 

Model. — In this Letter, we study a strictly ID model 
inspired by the Calogero-Sutherland (CS) model 12. 13|. 



mapping which states that the ground state wave 
function of the bosonic (B) and fcrmionic (F) states are 
related by *b(-R) = |*f(-R)|, where R = (ri,r 2 , . . . ,r„) 
are the one-particle coordinates. In case of bosons, the 



We consider the rotation-invariant ground state of n elec- 
trons on a ring of radius R. This yields a uniform electron 
density p = n/{2irR) = 1/(2 r s ), where r s = nR/n is the 
Seitz radius. This paradigm has been intensively stud- 
ied as a model for quantum rings, both experimentally 
[H - 21 j and theoretically 0, [22T - |29j because of their rich 
electronic, magnetic and optical properties, such as the 
Aharonov-Bohm effect and their potential application in 
quantum information theory. However, unlike the 1/x 2 
potential of the CS model, our electrons interact via the 
true Coulomb potential 1/x. 

Our work is inspired by three landmark papers. First, 
following Gell-Mann and Brueckner [3(|, we study the 
reduced (i.e. per electron) correlation energy, defined as 
the difference between the exact and Hartree-Fock (HF) 
reduced energies e c (r Sl n) = e(r a ,n) — e-nF{r s ,n), in the 
high-density (small-r s ) regime. We show that, despite 
using the true Coulomb operator, one can safely apply 
standard perturbation theory in this regime. Second, fol- 
lowing Wigner 31 1 , we study the correlation energy in 
the low-density regime, where the electrons crystallize 
to form a Wigner crystal. Combining the information 
obtained for these two limiting regimes, we propose a 
correlation functional that yields satisfactory estimates 
of the correlation energy for finite and infinite systems 
at high, intermediate and low densities. Third, following 
Ceperley and Alder [321 ], we explore the accuracy of the 
correlation functional by comparing its predictions with 
accurate diffusion Monte Carlo (DMC) calculations on 
finite and infinite systems. Wc also report an application 



2 



to a non-uniform system. Reduced energies and atomic 
units are used throughout. 

Hamiltonian. — The Hamiltonian (33| of n electrons on 
a ring of radius R is 



H 



1 



if 



2R 2 ^ dB 2 

i—1 1 



n 1 
E Z~> 



(1) 



»<7 



where #i is the angle of electron i around the ring cen- 
ter, and rij = R\J2 — 2 cos (6*j — 6j) is the interelectronic 
distance between electrons i and j. 

Hartree-Fock approximation. — The HF wave function 
^hf is a determinant of one-electron orbitals Xk(&) = 
e lke with orbital energy Kk = k 2 /(2R 2 ) : where k G Z 



if n is odd and k 



£ Z if n is even. These orbitals 



form a Vandermonde matrix [34j and, following the ap- 



proach of Mitas [35[ , one discovers the remarkable result 
*hf oc II" , r,,. where hj = 2i?sin[(^ - 6» J )/2] is a 
signed interelectronic distance. One sees that 'J'hf an- 
tisymmetric with respect to electron exchange and van- 
ishes whenever 9i = 6j . The resulting HF reduced energy 



enF(r s ,n) 



n 2 -l n 2 



A 4 - 1/n 2 



n 2 24r? 4 r 2fc - 1 

s s \fe=i 



3J, 

(2) 

where the first and second terms in (J2J represent the ki- 
netic and potential energies, respectively. Although the 
latter is finite for finite n, it cannot be partitioned into 
Coulomb and exchange parts, because each diverges. 

High-density expansion. — In the high-density 
{i.e. small r s ) regime, the kinetic energy is dominant and 
it is natural to define a zeroth-order Hamiltonian and 
a perturbation by = -1/(2R 2 ) £JU <9 2 /<90f and 

V = Y^i<j r iji an d the reduced energy expansion is 



s(r„n) = ^ + ^+e< 2 >(n) + 0(r a ). (3) 



The zeroth-order e^(n) = <*hf|^ (0) |*hf) and first- 
order eW(n) = (^hfI^I^hf) energies are found in ((5J. 
The second-order energy is 



A2) 
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K s K a Kb 



(4) 



where ty r A is a doubly-substituted determinant in which 
two electrons are promoted from {xa, Xb} to {xr, Xs}- If 
angular momentum is conserved, i.e. a + b = r + s, then 
K r + k s — K a — Kb — (r — a)(r — b)/R 2 , and the Slater- 
Condon rules reveal that, if a < b. 



k=r-b+l 



(5) 



Combining these yields 



(n) 



1 -1/2 

1 E E 

■ "- 1 b=a+l , 



^ (2-J- n , t )^ 6 , r 



(r — a)(r — 6) 



which can be evaluated to give, for example, 

e( 2 >(2) = 1 - 10/tt 2 , 

e (2) (3) = 16/9- 1436/(81tt 2 ), 

e (2) (4) = 109/45 - 244168/(10125tt 2 ). 



(7) 

(8) 
(9) 



Other values are shown in the second column of Table U 
and, in the r s — > limit, this gives all of the correlation 
energy. This many-electron system is one of the few for 
which one can obtain the exact closed-form correlation 
energy for any value of n. 

Low-density expansion. — In the low-density (r s > 2 
25 ) regime, the electrons form a Wiener crystal. Using 
strong-coupling perturbation theory [36j , the energy can 
be written 



r/°)(n) 



^j?W- 2 ), (10) 



where the first term represents the classical Coulomb en- 
ergy of the crystal and the second is the zero-point energy 
of the electrons vibrating around their equilibrium posi- 
tions. 

The Wigner crystal consists of n electrons separated 
by an angle 27r/n and is closely related to the one- 
dimensional Thomson problem. Thus, we have 



i<j 



(j - i)n 



(11) 



The second term in the expansion (JTUJ) is found by 
summing the frequencies of the normal modes obtained 
by diagonalization of the Hessian matrix. For electrons 
on a ring, the Hessian is a circulant matrix and its eigen- 
values and eigenvectors can be found in compact form, 
yielding 



r 3/2 



E 



2 - sin 2 



E 



4n 5/2 \ Z-j gin 3 tkv\ 
i=l \ k=l V n I 



sin 2 | % 



kn 



(12) 

Thermodynamic limit. — In the n — ¥ oo limit within 
the high-density regime, the kinetic energy 

£ ' » = tt 2 /24 (13) 

reduces to that of the ideal Fermi gas in ID @, H3 and 
the slow decay of the Coulomb operator causes 



= In Vn + (In 2 + 7 /2 - 3/4) + o(n°) 



(14) 



(where 7 is the Euler-Mascheroni constant [3J]) to grow 
logarithmically. 
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TABLE I. Reduced correlation energy (— e c (r a , n) in millihartree) for n electrons on a ring with various r s . Subscripts represent 
the statistical errors in the last digits. The DMC results for n — > oo are taken from Ref. Q. 
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TOT 


JJIVIO 


TOT 
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2 


13.212 


13.0o 


12.9 


12.77 


12.67 


12.15 


11.94 


11.250o 


10.922 


7.111 Oo 


6.733 6 


4.937 7o 


4.691 9 


3.122 02o 


2.998 22 


3 


18.484 


18.1 


18.0 


17.75 


17.64 


16.76o 


16.53 


15.346o 


14.997 


9.368 60 


8.972 2 


6.427 Oo 


6.168 8 


4.029 51o 


3.899 54 


4 


21.174 


20.5 2 


20.6 


20.24a 


20.15 


19.00i 


18.82 


17.320i 


17.003 


10.390 2i 


10.013 


7.084 5i 


6.838 2 


4.425 31i 


4.299 35 


5 


22.756 


22.3 2 


22.2 


21.662 


21.62 


20.33i 


20.15 


18.439i 


18.156 


10.946 3i 


10.590 2 


7.439 0i 


7.203 4 


4.635 56i 


4.514 64 


6 


23.775 


22.9 2 


23.2 


22.63 2 


22.57 


21.14i 


21.00 


19.137i 


18.887 


11.285 li 


10.947 3 


7.653 4i 


7.426 8 


4.761 89i 


4.645 07 


7 


24.476 


23.8 2 


23.8 


23.24 2 


23.21 


21.70i 


21.58 


19.607i 


19.384 


11.508 5i 


11.185 3 


7.794 81 


7.574 5 


4.844 00i 


4.730 64 


8 


24.981 


24.2 2 


24.3 23.69 3 


23.68 


22. Hi 


22.00 


19.940i 


19.739 


11.663 81 


11.352 9 


7.889 81 


7.677 7 


4.900 681 


4.790 10 


9 


25.360 


24.5 2 


24.7 24.04 2 


24.03 


22.39i 


22.31 


20.186i 


20.002 


11.776 5i 


11.475 8 


7.959 5i 


7.753 


4.941 36i 


4.833 27 


10 


25.651 


25. 2 4 


24.9 


24.25 4 


24.29 


22.62i 


22.55 


20.373i 


20.204 


11.856 81 


11.569 


8.013 4i 


7.809 8 


4.972 54i 


4.865 69 


CO 


27.416 




26.6 




25.90 




23.96 


21.444 1 2 


21.392 


12.317 74 2 


12.0913 


8.292 096 9 


8.120 5 5.132 504a 


5.039 13 



The limiting second-energy e^ 2 ' = — 7r 2 /360 can be 
found by converting the summations in ([6]) into integrals. 
In this way, one finds 



3 (r a ) = -0.027416 + 0(r s ), 



(15) 



and, in the dual thermodynamic/high-density limit, the 
exact correlation energy is therefore 27.4 millihartrees 
(mii'h) per electron. Using a quasi-lD model with a trans- 
verse harmonic potential, Casula et al. conclude that, in 
the high-density limit, the correlation energy vanishes 
Q. This strikingly different prediction stresses the im- 
portance of employing a realistic Coulomb operator. 
In the n — > 00 limit in the low-density regime, one finds 



(o) 1 r , ln ( 2 A) + 7 . , \ n(t , 
* > = m y/n H h o(n ), (16) 



n 



which has the same logarithmic divergence as euw, but 
with a different constant term. One can also show that 



(i) 



i- £ ^2Li 3 (l) - Li 3 (e«) - Li 3 (e-*)d0, (17) 



where Li 3 is the trilogarithm function |34| . We have not 
been able to find this integral in closed form, but it can 
be computed numerically with high precision, and yields 
r/ 1 ) = 0.359933, which is identical to the value found by 
Fogler in Ref. Q for an infinite ultrathin wire and a poten- 
tial of the form 1/ \J x 2 + /i 2 . This shows that, unlike the 
high-density limit where the details of the interelectronic 
potential are critically important, the correct low-density 
result can be obtained by using a modified Coulomb po- 
tential. Thus, in the dual thermodynamic/low-density 
region, we have 



e c (r s ) = -[ln(v / 27r)-3/4] r' 1 + 0.359933 r~ 3/2 + 0{r~ 2 ). 

(18) 



The same expansion can be derived for the infinite wire 
Q , confirming the equivalence of the electrons-on-a-ring 
and electrons-on-a-wirc models in the thermodynamic 
limit 
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Correlation functional. — We now use the information 
obtained in the high- and low-density limits to build a 
correlation functional for finite and infinite ID systems. 
We employ an interpolation between the high- and low- 
density limits inspired by the 'interaction-strength inter- 
polation" (ISI) expression of Seidl and coworkers [38( . We 
define 



TSI 



(r s ,n) 



t/W (n) [r? (1) (n)]< 



/3(n) In 



1/1 + a 2 (n)r 
^/l + a 2 {n)r s + (3{n)\ 



2e( 2 )(n)r 2 



l + /?(n) 



with 



o(n) = 2e (2) (n)7 ? (1) (n)/7 ? (0) (n) 2 
0(n) = e^(n)r]^(n) 2 /7]^(n) 3 - 



1 



, (19) 



(20) 
(21) 



Expression (fl9"|) reproduces the first term of the high- 
density expansion (fTS"|) and the first two terms of the 
low-density expansion (fTS]) . thus giving correct energies 
at r s = and vanishing at the correct rate for large r s . 

Discussion. — In TableHJ we report correlation energies 
for 2 < n < 10 and various r s values from the highl- 
and low-density regimes. Results obtained using the ISI 
method (|T9^) are compared with DMC calculations. Our 
DMC code follows the implementation of Ref. [39[ and 
the energy at zero time-step is obtained by linear ex- 
trapolation. The extrapolated standard error is obtained 
by assuming that the data are Gaussian distributed |4(J ■ 

The trial wave function is * T = *hf Yl7<j feiLi c k r ij) 
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and the coefficients have been optimized using the proce- 
dure described in Rcf. 
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For 0.1 < r s < 0.5, the ISI and DMC correlation ener- 
gies agree to within 0.2 mE^, which is remarkable given 
the simplicity of the functional. The deviations increase 
to 0.3-0.4 nxEh for r s = 1 and r s = 5, but then decrease 
again to 0.1-0.2 mE^ for r s = 10 and r s = 20. Overall, 
the ISI correlation functional gives reasonable estimates 
of the correlation energy for 2 < n < 10. 

The functional can also be used to estimate the corre- 
lation energy in infinite systems and, in Table [H its pre- 
dictions are compared with the DMC calculations of Lee 
and Drummond Q for an infinitely thin wire. We find 
that e* SI underestimates the energies by 0.05, 0.23, 0.17 
and 0.09 mE^ for r s = 1, 5, 10 and 20, respectively, indi- 
cating that Eq. (fl9|) provides accurate correlation energy 
estimates in both finite and infinite uniform systems. 

Moreover, the expression (fT9|) can also be used as a 
correlation functional for non-uniform systems. For ex- 
ample, let us consider a ID two-electron quantum dot 
(QD), i.e. two electrons interacting via a bare Coulomb 
potential and confined by a harmonic potential of force 
constant a; 2 . Following Ref. [2^|, we find that, for 
ui = 1/2, the exact reduced energy is e QD = 3/4, the 
exact wave function is ip® (xi,X2) = (x\ — £2)(1 + 
\ x i — a?2 1/2) exp[— {x\ + x|)/4] and the exact density is 
p QD (x) = {7+x(2+.T)[6+a;(2+a;)]}exp(-a; 2 /2)/(10v / 2^ : ), 
where xi and xi are the Cartesian coordinates of the 
electrons. Using a Gaussian expansion, we obtain e§p = 
0.7591, which yields e Q D = -0.0092. The expression 1(15]) 
with n = 2 yields 



QD [p QD l - 



1 



p QD {x)ef I (p QD ,2)dx = -0.0106, 



(22) 

which underestimates the reduced correlation energy by 
only 1.4 m£"h. We note that if one use the functional (fT9)) 
with n — > oo, the resulting correlation energy is found to 
be wrong by a factor two. 
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